Pressure Induced Complexity in a Lithium Monolayer 
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Light alkali metals have usually been considered as simple metals due to their monovalency and 
high conductivity. In these metals ionic pseudopotentials are weak and the nearly free electron model 
(NFE) becomes quite accurate at normal conditions. However, very recent experiments have shown 
that at high pressures their electronic properties deviate radically from the NFE model and even 
become unexpected good superconductors. In this work we present ah initio calculations to analyze 
the deviation from simplicity in a lithium monolayer (ML) when pressure is applied. We have seen 
that as a result of the increasing non-local character of the atomic pseudopotential with increasing 
pressure, the surprising half filling Hubbard-type nesting observed in the Fermi line can explain 
the interesting complex behavior in lithium ML, induced by its correlated structural, electronic and 
even magnetic properties. 

PACS numbers: 71.18.+y, 71.15.Mb, 68.35.-p, 73.20.At 



I. INTRODUCTION 

Light alkalies are usually considered as simple metals, 
because as a consequence of the weak interaction 
between valence electrons and ionic cores, under normal 
conditions they crystallize in simple high symmetric 
structures. Therefore, the nearly free electron model 
(NFE) has been considered a good approximation 
to describe their electronic properties The only 

exception is hydrogen which, even when solidifies, 
forms very stable diatomic molecules and presents an 
insulating character. However, recent both theoretical 
and experimental results have changed our perception of 
the alkalies and have shown that the NFE model breaks 
when high pressures are applied. Neaton et al. reported 
ah initio calculations for bulk lithium P| and predicted 
that pressure could induce structural transitions to less 
symmetric, lower coordinated structures, associated with 
electronic localizations. These theoretical predictions 
have been confirmed by experiments finding that 
around 40 GPa a phase transition to a complex structure 
(c/16) with 16 atoms per unit cell occurs. It is important 
to notice that pressure induced transitions from simple 
to more complex structures are not singular to lithium 
but have also been predicted in heavier alkalies SpL On 
the other hand, other very recent experiments 0,0, |1| 
have also shown that at this pressure range lithium 
presents a superconducting transition with Tc 20K, 
becoming the highest transition temperature between 
simple elements |9|. It is noteworthy that experiments 
looking for superconductivity in lithium under ambient 
pressure have failed This even rises the interest 

to characterice the physical properties of compressed 
lithium close to the observed electronic and structural 
transitions. 

In this article we present an ah initio theoretical analysis 
of the pressure-induced structural, electronic and even 
magnetic properties of a lithium ML. The develop- 



ment of new techniques for the atomic manipulation 
[Til allows the growth of atomic monolayers on 
inert substratum, semiconductors or noble gases. These 
new possibilities rise the interest to analyze physical 
properties of low dimensional systems under different 
conditions. On the other hand, the simplicity of the 
atomic configuration associated to the ML facilitates 
the performing of a more detailed analysis, and extend- 
ing our conclusions to the bulk will also provide an 
interesting perspective to understand the physical origin 
of the experimentally observed features of compressed 
lithium. In Section ^ we describe the theoretical and 
computational background of this work. Results and 
Conclusions will be presented in Sections 11111 and IIVI 
respectively. 



II. THEORETICAL AND COMPUTATIONAL 
BACKGROUND 

In order to describe structural and electronic properties 
of a lithium ML over a wide range of densities we have 
used a plane-wave implementation of the Density Func- 
tional Theory (DFT) fll. [11 w ithin the Local Density 
Approximation (LDA) [lllll 13 

for non magnetic 

calculations and Local Spin Density Approximation 
(LSDA) to analyze the magnetic properties of the 
ML. All calculations have been performed with the 
Vienna Ab initio Simulation Package (VASP)[l3, H^. 
Effects of compression were simulated by reducing the 
lattice parameter. As the atomic spacing is reduced 
a significant core overlap is observed, so that we have 
required methods which take into account core electrons. 
For fully treating them we make use of the Projector 
Augmented Wave (PAW) method HJ. A Monkhorst- 
Pack 22] mesh of 20 x 20 x 1 has been chosen for the 
sampling of the Brillouin zone. The description of the 
ML is implemented by a supercell which, in order to 



2 



minimice interactions between MLs, contains ten layers 
of vacuum between MLs. The ions are relaxed to their 
equilibrium positions by calculating the forces acting 
on them using the Hellmann-Feynmann theorem. In 
all the calculations we have used an energy cutoff of 
-Bcut ^ 815eV. 



at r decreases so much that the Pz band, antisymmetric 
in the direction perpendicular to the ML, starts to be 
occupied and opens a new conduction channel. It is 
also interesting to note that the van Hove singularities 
lay close to the Fermi energy which, therefore, rises the 
DOS. All these interesting characteristics can only be 
explained in terms of the non-local character of the ionic 
pseudopotential ,23J. 



III. RESULTS 
A. Structural and Electronic Analysis. 

Following the analysis in Ref. |23] , as a first step we will 
just consider three simple structures of a lithium ML: 
hexagonal (hex), square {sq) and honeycomb (he). Fig. 
^shows the calculated enthalpies for these structures. As 
it is a common property of the alkalies, lithium ML at 
equilibrium (r^ — 3.02 a.u.) |2^ also favors the compact 
hexagonal structure (hex). However, as we can observe 
in Fig^ at < 2.25 a.u. the more open square structure 
(sq) becomes favored over the hexagonal one. 
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FIG. 1: Enthalpy H = E + a A, a being the surface tension: 
a — —dE/dA, as a function of the density parameter, r^, 
for three simple structures: hexagonal (hex), square (sq) and 
honeycomb (he). All enthalpies are referenced to the square 
structure, AH = H — Hsq. 

In Fig. |3we plot the band structure and density of states 
(DOS) of a lithium ML corresponding to the hexagonal 
structure at equilibrium and the square structure with 
Ts = 2.15 a.u., just after the structural transition. At 
equilibrium (Fig. 2a) only one band is occupied, which 
is characterized by a parabolic-type dispersion of s 
character, as corresponds to a good NFE-like system. 
For occupied energies the DOS is almost constant, which 
reflects the quasi two-dimensional character of the ML. 
However, at higher pressures (Fig. 2b) the first band 
flattens associated to the increasing of the band gap at 
the zone boundary, indicating a clear electronic localiza- 
tion. The energy difference between the first two bands 




FIG. 2: Band structure and density of states (DOS) for (a) 
a hexagonal lithium ML at equilibrium (rs = 3.02 a.u.) and 
(b) a square structure with rs = 2.15 a.u. The Fermi energy 
is represented by the dotted line. 

Another essential feature determining electronic prop- 
erties of the ML is the geometry associated to the 
Fermi line, which is displayed in Fig. U] for the same 
densities and structures as above. At equilibrium just 
one Fermi line exists, associated to the s band, which is 
almost circular as corresponds to a NFE approximation. 
However, as mentioned above, this simple picture breaks 
at higher densities. At rg = 2.15 a.u. there are two 
Fermi lines associated to the two occupied s and Pz 
bands and, surprisingly, the Fermi line associated to 
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the s band becomes a perfect square as corresponds to 
a half filled Hubbard- type model. As it is well known, 
this model becomes appropriate to describe electronic 
properties of systems with localized electronic states, 
just being the opposite case to the NFE approximation 
valid at normal pressures. The perfect nesting observed 
at rg = 2.15 a.u. strongly couples electronic states close 
to the Fermi line along the TM direction, and might 
induce structural Peierls transitions. 
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FIG. 3: Fermi hues of a) a hexagonal lithium ML at equilib- 
rimn and b) a square ML with = 2.15 a.u. The solid line 
corresponds to the Fermi line of the s band and the dotted 
line is associated to the Pz band, which becomes occupied at 
high densities. 



B. Peierls Transition 

As a consequence of the nesting in the Fermi line at 
rg ~ 2.15 a.u., any perturbation which breaks the 
electronic degeneration might be favored. Although 
different symmetry breaking perturbations might be 
proposed, in this section we have included the analysis 
of the two simplest diatomic square structures: dsqi and 
dsq2. In both structures we choose a diatomic square 
unit cell and allow the central atom to move along 
two highly symmetric directions: the nearest neigh- 
bor (dsqi) and second nearest neighbor {dsq2) directions. 

In Fig. ^we present the enthalpies as a function of the 
density for the five structures considered: he, hex, sq, 
dsqi and dsq2. It is interesting to note that at ~ 2.15 
a.u, just the density where the nesting in the Fermi line 
of the monoatomic square structure (sq) is found, the 
new dsq2 structure proposed in this work is favored over 
the others. Therefore, the perturbation caused by the 
nesting in the Fermi line consists in a relative movement 
of the atoms in the second nearest neighbor direction, 
producing a charge density wave in this direction. 
However, at much higher densities (r^ — 1.35 a.u.) the 
sq structure is favored again over the diatomic structures. 

Fig. |S1 shows the band structure and DOS for the dsq2 
structure at = 2.0 a.u., which confirms the clear 
deviation from a NFE behavior, already manifested for 
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FIG. 4: Entalphy as a function of the density parameter, 
rs, for all the structures considered: he, hex, sq, dsqi and 
dsq2. All the enthalpies are referenced to the square enthalpy, 
AH = H — Hsq- Although the sq structure becomes favored 
at Vs = 2.25 a.u., the proposed new diatomic square structure 
{dsq2), where the atom at the center is allowed to move along 
the direction of the second nearest neighbor, is preferred at 
higher densities (1.3 a.u.< <2.15 a.u.). 

the monoatomic structures of a lithium ML. The s band 
amazingly flattens in the XM direction, which indicates 
a strong electronic localization in this direction, although 
the ML remains metallic. In Fig. we plot the LDA 
valence charge density of the dsq2 structure at — 2.0 
a.u. The charge density accumulates between the ions 
forming a clear one-dimensional conducting channel, 
which also originates a very sharp singularity in the 
electronic DOS close to the Fermi energy (Fig. |3J). 
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FIG. 5: Band structure and DOS of a lithium ML with dsq2 
structure at = 2.0 a.u. The s band becomes extraordinary 
flat in the XM direction. 

In order to characterize the phase transition involving 
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FIG. 6: LDA valence charge density of the proposed dsq2 
lithium ML at Vs = 2.0 a.u. The ions with core electrons 
are represented by circles with dots inside, while the valence 
charge density maximum (black areas) accumulates between 
them forming a clear one-dimensional conducting channel. 



the dsq2 structure we have analyzed the displacement of 
the central atom along the minimum energy direction 
as a function of the density. As it is shown in Fig. 
H and expected from our analysis above, the atomic 
displacement is zero for rs > 2.15 a.u. However, at 
Ts = 2.15 a.u. the high anharmonic potential of the 
central atom induces the presence of a soft phonon mode 
which moves the atom at the center. As it is displayed 
in Fig. [7| around — 2.15 a.u. the displacement of 
the central atom grows continuously with increasing 
density which indicates the second-order character of 
the sq — > dsq2 transition. However, if we keep rising the 
density, we find that at ~ 1.3 a.u. the displacement 
of the central atom vanishes discontinuously and the 
ML goes back to the monoatomic sq structure. We have 
seen that this transition is explained by the presence 
of an absolute minimum energy at the center while 
relative minima at finite displacements from the center 
still remain, characterizing the dsq2 sq transition as 
first-order. 



C. Magnetic Instability 

In a low density electron gas it is well known that a ferro- 
magnetic instability can be induced as the kinetic energy 
penalty in the ferromagnetic state is compensated by the 
gain associated to the electronic exchange l^^. So that, 
in principle, we would not expect such a transition to oc- 
cur for high densities. However, as we can see in Figs. [21 
and [SI the second band starts to be occupied as density 
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FIG. 7: Displacement of the central atom in the dsq2 structure 
(solid dots) along the second nearest neighbor direction as a 
function of the density. The open squares correspond to the 
sq structure, where the atom remains at the center. 



is increased. Therefore, the effective low electronic den- 
sity in the Pz band, in conjunction with the high DOS 
in the proximities of the Fermi energy and the nesting 
observed in the Fermi line, could induce a ferromagnetic 
instability in a lithium ML . Fig. [Sjdisplays ab initio cal- 
culations within the LSDA approximation which predict 
that a ferromagnetic instability can occur in a narrow 
density range: 2.2 a.u.> > 1.95 a.u. It is interesting 
to notice that the maximum value of the net magnetic 
moment per conduction electron, 0.14/iBj is observed 
at Ts =2.1 a.u, just at the density where the sq dsq2 
Peierls structural transition takes place. This magnetic 
instability in a lithium ML, which has been shown to be 
correlated with electronic and structural transitions pre- 
sented above, also proves the complexity that pressure 
induces in these, otherwise considered, simple systems. 



IV. CONCLUSIONS 

In summary, according to the ab initio calculations 
presented in this work, although the NFE model cor- 
rectly describes the properties of the lithium ML at 
normal conditions, it loses its validity as the electronic 
density increases. We have analyzed pressure induced 
complexity in structural, electronic and even magnetic 
properties of the ML, finding an interesting correlation 
between them. As density goes up, a flattening of the s 
band is observed, indicating an electronic localization. 
The corresponding nesting in the Fermi line and van 
Hove singularities in the proximities of the Fermi energy 
observed at Ts ~ 2.1 a.u. induce a Peierls second-order 
structural transition {sq dsq2). Beyond these features, 
at around the same density, the non-local character of 
the ionic pseudopotential allows to start filling the pz 
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FIG. 8: Magnetic moment per atom as a function of the 
density, Vs, for a lithium ML with sq (squares) and dsq2 
(dots) structures, which minimize the energy at the pressure 
range considered. Lithium ML presents a stable ferromag- 
netic ground state in the proximities of the Peierls transition, 
with a maximum net magnetic moment per atom of ~ 0.14/iB, 
located close to the density where the structural transition 
(sq dsq2) has been predicted, = 2.1 a.u. 



band, and its low effective initial occupation causes a 
ferromagnetic Bloch instability. Besides of the intrinsic 
interest of studying the physical properties of systems 
with reduced dimensionality, the simplicity of the ML 
allows us to make a more detailed analysis and by 
extending these conclusions to bulk systems will provide 
us a very useful perspective to understand the physical 
origin of the complexity observed in compressed bulk 
lithium [26l| . 
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